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Abstract 

One of the crucial differences between mathematical models of classi- 
cal and quantum mechanics is the use of the tensor product of the state 
spaces of subsystems as the state space of the corresponding composite 
system. (To describe an ensemble of classical composite systems one uses 
random variables taking values in the Cartesian product of the state spaces 
of subsystems.) We show that, nevertheless, it is possible to establish a 
natural correspondence between the classical and quantum probabilistic 
descriptions of composite systems. Quantum averages for composite sys- 
tems (including entangled) can be represented as averages with respect to 
classical random fields. It is essentially what Albert Einstein was dreamed 
of. Quantum mechanics is represented as classical statistical mechanics 
with infinite-dimensional phase space. While the mathematical construc- 
tion is completely rigorous, its physical interpretation is a complicated 
problem (which will not be discussed in this paper). 



1 Introduction 

Nowadays it is commonly accepted that the use of the tensor product of 
the state spaces of subsystems as the state space of the corresponding 
composite system is one of the main distinguishing features of QM. It 
is especially important in quantum information theory where systems in 
entangled states play an fundamental role0 In this paper we do not discuss 
extremely complicated problems related to interpretations of quantum 
mechanics. We proceed in the framework of mathematical physics. 

There are known two models [S] for computations of averages for en- 
sembles of composite systems: a) classical probability model (due to Kol- 
mogorov [9]) based on integrals, b) quantum probability model (due to 
von Neumann [10], see, e.g., |12) . |13) for the modern treatment of 

the problem) based on traces of self-adjoint operators. We show that, in 

^ While it is sufficiently well studied mathematically, entanglement is still quite mysterious 
physically. Its widely used interpretation as the evidence of "nonlocal correlations" can not be 
considered as completely satisfactory, see, e.g., [T]- [5] and especially [6]- [8] for recent debates. 
Therefore clarification of the structure of its mathematical description may have important 
consequences. This paper is a step in this direction. 



spite of a rather common opinion, quantum correlations for observables 
on subsystems of a composite systems can be represented as correlations 
with respect to a classical (Gaussian) random fields. Moreover, dynamics 
of quantum correlations induced by Schrodinger's equation can be re- 
duced to dynamics of correlations for a classical "prequantum" stochastic 
processes. 

For non-composite systems theory has been developed in a series of au- 
thor's papers |14) - [2U| . It is known under the name Prequantum Classical 
Statistical Field Theory (PCSFT). In this theory ensembles of quantum 
particles are represented by classical random fields, probability distribu- 
tions on a Hilbert space. We remark that appealing to classical random 
fields is rather common in various attempts to create a kind of classical 
statistical mechanics which reproduces predictions of QM. We can men- 
tion stochastic electrodynamics, e.g., [21) . |25) . or the semiclassical model, 
e.g., |26|- [28] . Bohmian mechanics also contains a kind of classical field, 
the pilot wave. However, in this model randomness is coupled to particles 
and not to fields. The same can be said about Nelson's stochastic QM 
[29| and its generalization due to Mark Davidson |30| . [31| as well as the 
recent prequantum model of 't Hooft [32| . [33| . see also Thomas Elze [35]. 
Physically the present paper belongs to the domain of "quantum mechan- 
ics as emergent phenomenon", cf. mentioned papers of 't Hooft and Elze 
as well as models which were recently created Rusov et al. '36] and Kisil 



However, as was already pointed, we prefer not to go in the debate 
on a physical meaning of the proposed mathematical construction. Our 
aim was to unify two mathematical descriptions of averages, classical and 
quantum. This aim was approached via representation of quantum cor- 
relations by Gaussian integrals over the Hilbert space Hi x H2, where 
Hi.i — 1,2, are the (Hilbert) state spaces of the subsystems. On the 
other hand, we could not totally escape the interpretation problem. The 
main message from our mathematical construction is that a quantum pure 
state of a composite system should be considered not as a "state vector" 
belonging to the tensor product Hi ® H2, but as non-diagonal block of 
an operator acting in the Cartesian product Hi x H2 see (|35|l . This 
operator, say D, is the covariance operator of the prequantum classical 
random field. We remind that consideration of a density operator as the 
covariance operator of the corresponding prequantum random process is 
the crucial point of PCSFT, see |14) - |20| . In this paper we extend this 
approach to composite systems. 

We point out to another approach providing a possibility to represent 
quantum averages by operating only with classical probability distribu- 



tions, namely, quantum tomographic approach, see Manko et al. [38] 



Results of these paper were shortly announced in [43], where they 
were presented on the physical level of rigorousness. Unfortunately, such 
a formal presentation induced (to be totally honest) a mistake - unfortu- 
nately, in the basic equality (|22[l coupling the quantum correlation with 
the nontrivial term of the classical field correlation. In this paper we 
proceed in the rigorous mathematical framework, the basic formula and 
its consequences were corrected. Formulas of the present theory are more 
complicated than of one given in ,43j . but they are mathematically correct. 
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2 Statistical models of classical and quan- 
tum mechanics 



Everywhere in this papers Hilbert spaces are separable. Let Hj,j = 1,2, 
be (real or complex) Hilbert spaces. We denote the space of bounded 
linear operators from Hi to H2 by the symbol £{Hi, H2)- Let _ff be a (real 
or complex) Hilbert space. We denote the space of self-adjoint bounded 
operators in H by the symbol Ca[H). 

2.1 Classical model 

a). States are represented by points of some set M (state space). 

b) . Physical variables are represented by functions / : M — >■ R be- 
longing to some functional space 

c) . Statistical states are represented by probability measures on M 
belonging to some class S{M). 

d) . The average of a physical variable (which is represented by a func- 
tion / e V{M)) with respect to a statistical state (which is represented 
by a probability measure p G S{M)) is given by 

< / >p= / fmp{<p). (1) 

JM 

By using the language of probability theory we can say that there is given 
a random vector <^(a;), where w is a random parameter, taking values in 
M. Then {f}^ — Ef[(j>{ui)) — {f)p. Here and everywhere below E denotes 
classical mathematical expectation (average). 

If the state space M is a space of functions, e.g., M = L2(R^), then 
M- valued random vectors are called random fields. For each cj, <^(a;) is a 
function of a; £ R'^ : 4>{x,uj). 

e) . If systems Si,i = 1, 2, fc, have state spaces Mi, respectively, then 
the composite system S — {Si)i^i has state space M = Mi x . . . x Mk, 
the Cartesian product of the state spaces Mi. Ensembles of S-systems 
are described by random vectors in M : (t){ijj) — {<j)i{u)) , . . . , (j>k{oj)) (or 
equivalently by probability measures on M.) A trivial, but important, 
remark is that in general components of ^(a;) are not independent. There 
are nontrivial correlations between them. The best way to describe these 
correlations is to use the covariance operator (it will be defined little bit 
later). 

A classical statistical model is a pair Af = (S(M), V{M)). 

2.2 Quantum case 

Let H he a complex Hilbert space. 

a). States (pure) are represented by classes of normalized vectors of 
H with respect to the equivalence relation: Tpi = e'®i/'2. 

a). Physical observables are represented by operators A : H ^ H 
of the class Cb{H). (To simplify considerations, we shall consider only 
quantum observables represented by bounded operators.) 

^The choice of a concrete functional space V{M) depends on various physical and mathe- 
matical factors. In classical mechanics for systems with the finite number of degrees of freedom 
M is chosen as the phase space R^"; V(M) is the space of smooth functions. 
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b). Statistical states are represented by density operators. The class 
of such operators is denoted by 'D{H). 

d) . Average of a physical observable (which is represented by the op- 
erator A G £s{H)) with respect to a statistical state (which is represented 
by the density operator p £ D{H)) is given by von Neumann's formula 

< A >p= Tr pA (2) 

e) . If quantum systems Si,i = 1,2, ...,k, have the state spaces Hi, 
respectively, then the system S — {Si)i^i has the state space Hi0. . .(E)Hk, 
the tensor product of state spaces Hi. 

The quantum statistical model is the pair A^quant = {'D{H), Cs{Hy). 

At the first sight the gap between classical statistical mechanics and 
quantum mechanics is huge [10| . Impossibility to reduce quantum av- 
erages to classical averages is the main source of the great ideological 
difference between classical and quantum probabilistic descriptions. 



3 Gaussian measures on real and complex 
Hilbert spaces 

3.1 Real case 

Let W be a real Hilbert space. Let A £ Ca{W). 

We start with derivation of the basic mathematical formula which was 
used in [14| - |20| . We will calculate the Gaussian integral of the quadratic 
form 

fA{4>) = {Adf,4>). (3) 

Consider a cr-additive Gaussian measure p on the cr-field of Borel subsets 
of W. This measure is determined by its covariance operator B : W ^ W 
and mean value m £ W. For example, B and m determine the Fourier 
transform of p 

Jw 

In what follows we restrict our considerations to Gaussian measures with 
zero mean value: (m,y) = J^{y,^)dp{ip) = for any y £ W. Sometimes 
there will be used the symbol ps to denote the Gaussian measure with the 
covariance operator B and m — 0. We recall that the covariance operator 
B is defined by its bilinear form 

{Byi,y2) = J{yi,<j)){y2,<l))dp{(l)),yi,y2 £W, (4) 

and it has the following properties: a) B > 0, i.e., {By,y) > 0,y £ W; 
b) -B is a self-adjoint operator, B £ CsiW); c) B is a trace-class operator 
and Tr _B = \ \<f)\f dp{<f)) . It is dispersion of the probability p. Thus for 
Gaussian probability we have cr^(p) = Tr B. We remark that the list of 
properties of the covariation operator of a Gaussian measure differs from 
the list of properties of a von Neumann density operator only by one 
condition: Tr p = 1, for a density operator p. Thus, for any covariance 
operator B, its scaling B/TrB can be considered as a density operator. 
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By using Q we can easily find the Gaussian integral of the quadratic 
form fA{(t>) defined by ([3}: 

fA{<P)dpB{cb) = f {Acb,4>)dpB{<l>) 

Jw 

oo p oo 

= y'(^e!.ej) / {ei,<l>){ej,(l>)dpB{tp) ^ V (yle^, ej)(Bei, ej), 

where {e^} is some orthonormal basis in W. Thus fA{4')dpB{4') = 
Tr BA. 



3.2 Complex case 

Let Q and P be two copies of a real Hilbert space. Let us consider their 
Cartesian product H — Q x P, "phase space," endowed with the symplectic 

operator J = ^ ^ ^ ) ' class of Gaussian measures (with 

zero mean value) which are invariant with respect to the action of the 
operator J; denote this class S{H). It is easy to show that p G S{H) if 
and only if its covariance operator commutes with the symplectic operator, 

mi- 

As always, we consider complexification of H (which will be denoted 
by the same symbol) , H = Q(B iP- The complex scalar product is denoted 
by the symbol {■,•)• We will also use the operation of complex congugation 
* in complex Hilbert space H, for (j> = 4'i + i<l>2,(t>i G Q,02 £ P, we set 
*{4)) = (j> = (pi — i4>2- We will use the following trivial fact; 

{u,v)^{v,u). (5) 

We introduce the complex covariance operator of a measure p on the 
complex Hilbert space H 

{Dyi,y2)= / {yi,4>){(t>,y2)dp{4>). 

J H 

We also consider the complex Fourier transform of p 

p{y)=l exp{i{{y,(l>) + {4>,y))}dp{4>). (6) 

J H 

Any J-invariant Gaussian measure on H is determined by its complex 
Fourier transform. [44|: p(y) = exp{ — {Dy,y)}. 

We remark that J-invariance is a strong constraint on the class of 
Gaussian measures under consideration. Consider a measure p on the 
Cartesian product H = Q x P. Its real covariance operator has the block 
struct ur«[f| 

-B = I 1 , where B^i = i3ii,_B|2 = 522,-8^2 = B21. Con- 



-821 B22 

sider also its complex covariance operator D. It can be realized as acting 
in the Cartesian product of two real Hilbert spaces and in such a represen- 

L C 

tation it also has the block structure -Drcai ~ I ^ r 1 . It was shown 



-C L 

in 144) that L = Bn + B22 and C = -612 — -621- It also was shown that 



^Little bit later we will use the block structure of the complex covariance operator for a 
measure defined on the Cartesian product of two complex Hilbert spaces. The reader should 
be careful and not mix these two totally different block structures! 
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if measure is symplectically invariant then Bn = B22, B21 = —-812- Thus 
in the latter case the complex and real covariance operators are coupled 
in a simple way: Drcai = 2_B. 

Lemma 1. For any measure p G S{H) the following representation 
takes place 

/ (?i,<^>(77i,0>(0,6>(0,'72)dp(0) = {D^um){Driui2) + {Dii,i2){Drium) 

J H 

(7) 

To prove this formula, one should differentiate the Fourier transform 
^ four times. 

Let H\ and H2 be two complex Hilbert spaces and let D21 G C{Hi, H2), -D12 G 
C{H2, Hi). Then D21 ® -D12 G C{Hi (gi H2, H2 i/i). Let us consider the 
permutation operator a : H2 ® H\ Hi (g) H2, (j{4>2 ® 0i) = 4>i® 4>2- We 
remark that cr G C{H2 ® Hi, Hi (3 -^2)- 

Let p be a measure on the Cartesian product Hi x H2 of two Hilbert 
spaces. Then its covariance operator has the block structure 

where Da : Hi ^ Hi and Dij : Hj — > Hi. The operator is self-adjoint. 
Hence D*^ = Da, and _Di2 ~ D21. 

Let // be a complex Hilbert space and let A G C{H,H). We con- 
sider its quadratic form (which will play an important role in our further 
considerations) 

(l>-^ fA{4>) = {MA)- 

We make a trivial, but ideologically important remark: /a : -ff — > -ff, is a 
"usual function" which is defined point wise. 

In the same way as in the real case we prove the equality 

/ ,fA(0)dpo(0) = TrDA (9) 

J H 



Theorem 1. Let p G S{Hi x H2) with the (complex) covariance 
operator D and let Ai G C{Hi, Hi), i = 1,2. Then 

[ .fA^ {4>i).fA2 {4>2)dp{<j}) = TrDiili Tr7:i22 A2 + TrDi2A2D2iAi 
J H1XH2 

(10) 

This theorem is a consequence of the following general result: 
Lemma 2. Let p G S{H) with the (complex) covariance operator D 
and let Ai G £{H, H),i = 1,2. Then 

fA^{(l>).fAA<P)dp{<l>),= TrDAiTiDA2 + TiDA2DAi. (11) 

H 

Proof. By Lemma 1 the integral can be represented as 

/ = Gil, ejj (^26^2, 6^2) 

x[(-Deji,ei2)(L»ej2,eii) + (Dej-j , e.^ ) (De^^ , e,^)] ^ h + I2, 
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where {e,} is an orthonormal basis in H. Here 



»li2 jl 32 

i 

h = ^^{lieiieji)(eji,Deii)^{l2ei2ej2){ej2,£)ei2) 

i2 J2 

= ^(Aie,^ , Dci^ ) (lae,, , Da^ ) = TtDAiTtDa^ ■ 

ili2 

Proposition 1. Let conditions of Theorem 1 hold. Then 

/ai (0i)/a2 {(j>2)dp{4>) = Tr(Dii ig) D22 + ffCIJai ® -D12)) Ai ® A2. 

H1XH2 

(12) 

Proof. It is sufficient to prove tliat 

Tro-(D2i » Di2)Ai » I2 = T1D12A2D21A1. 

We have 

Tr(7(£)2i ® -Di2)Ai (g) A2 = ^(a(D2i ® -Di2)li ® Ae; ® /j, a ® /,) 
= Y,{Di2A2fj(giD2iAiei,ei ® /,) = ^(Dia^a/j, e.)(D2i^iei, /j). 

ij i3 

On the other hand, TrDi2A2D2iAi 
= ^{D-,2A2D2iMeuei)^{D2iMei,AlD2iei)^{D2iMeiJj){fj,AlD^ 

i i ij 

4 Vectors, operators, traces 

4.1 Vector and operator realizations of the tensor 
product 

In quantum theory a pure state of a composite system is represented by a 
normalized vector belonging to the tensor product Hi®H2. On the other 
hand, in functional analysis it is common to use elements of H\ (g) H2 
as operators. The standard construction provides the realization of any 
vector e (g) //2 by a linear operator from H2 — > i^i, where H2 is 
the space dual to H2. For a vector Z = u ® v, one puts Zx* = x*(v)u 
for X* € H2- This correspondence Z ^ Z is extended to isomorphism 
of Hi (g H2 with the space of Hilbert-Smidt operators 'HS{H2,Hi) or 
(equivalently) with the space of anti-linear ffS'-operators from H2 to Hi. 
The main compication in coming considerations is related to our need to 
represent vectors from the tesor product by linear operators from H2 to 
Hi. We present a new construction whicli seems to be unknown, so it may 
be interesting even from the purely mathematical viewpoint. We remind 
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the definition of tlie HS-noTva corresponding to tiie trace scalar product 
in 'HS{H2, Hi). Take an arbitrary orthonormal basis {fk} in H2- Then 

(Vi, T2) = ^{Ll/fc, L2h) = TrQ'Ll, ||L||^ = TrL*L 
fe 

(index 2 is typically used for the Hilbert-Schmidt norm). 

Let {e^} and be two orthonormal bases in H\ and H2, respectively. 
Then 

= Y^iJ^je^® fj,^lj,j ec, 

and ll^lP ~ IV'ijP- We remark that, for an orthonormal basis, say 
{fj}, the system of complex conjugate vectors {fj} is also an orthnormal 
basis. Really, by (O we get (f^J,) = (/,, fj) = Sij. 
We set, for </!> G H2, 

*'/' = EV'»i('^2j.>- (13) 

We emphasize that, for the vector 'i', the expansion with respect to the 
basis {ei(g)/j} was used. In cotrast to this, in the expansion of the operator 
^I" the basis {e; ® fj} was used. If the basis is real than the definition is 
essentially simplified: 

*0 = E'/''^-<'^2,/,)e.. (14) 

Take now a factorizable vector ^ = it ® u. It defines the rank one 
operator 

^(p = {(p, v)u — {v, (f))u. 
To play with the definition, take = cu®v = (cu) Cg) u = u {cv) , where 
c £ C. On the one hand, we have 'i'tf) = c{(f),v)u; on the other hand, 
'^'(f) — {(l},cv)u = {(j),cv}u — c{(f),v)u. Thus our definition is consistent 
with scaling by a complex constant. Consider now very special, but at 
the same time very important case: Hi = L2(R"'),i = 1,2. Here our 
definition gives the following representation (a special case of (|14|) 1: 

*0(a;) = j <^{x,y)ci>{y)dy (15) 

We proceed formally and use the real basis {ex®fy}, where exit) = 5{t—x) 
and fy{s) = S{s ~ y). Finally, we remind a result of theory of integral 
operators. The operator given by (|15|) is of the HS-type in L2-spaces 
if and only if its kernel ^{x,y) is square integrable. Thus the condition 
* e 1/2 (R"^ X R"2) for a function 5' is equivalent to the iyS-condition 
for the operator 5'. In the abstract form this fact will be formulated in 
coming lemma. 

Lemma 3. Each vector ^ G Hi ® H2 determines (uniquely) opera- 
tor ^ G C{H2,H{) and, moreover, ^ G HS{Hi,H2) and ||$||2 = ||*|| 
(the norm of a vector coicides with the Hilbert-Schmidt norm of the cor- 
responding operator) 

Proof. We start with proof of corectness of definition 1)13^ . i.e., it does 
not depend on the choice of orthonormal bases. Let {e'j} and {/j} be two 
orthonormal bases in Hi and H2, respectively, which are in general differ- 
ent from bases {cj} and {fj}- Consider the expansion ^ = "^-j ^l)[je'i® f'j. 
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We remark that /j = J^kifj' fk)fk = J2k(fk^f'j)fk, which can also be 
obtained by complex conjugation from the expansion /j = '}2k(f-i'fk)fk- 
In this basis the vector ^ defines the operator 

$> = ^V^:,(0,7J)e:. (16) 

ij 

We have 

> (/- , > (0, E </i ' ) ^ ) E ^j') 

nm kp i j 

We remark 

Thus ^'(/> = ''/'i™.(</') It coincides with (fT5)) . We now apply 

the Cauchy-Bunyakovsky inequality to the expression (|13p and obtain: 

Thus this operator belongs to the space C{H2, Hi). Finally, we show that it 
belongs even to the space HS{H2, Hi) : = Tr$*$ = Ek(*fk, *fk) = 

Ek(Eii V'iikBii.Eii -^iskeis) = Eki^ikV'ik- 

Moreover, it is easy to show that any operator L G 'HS{H2, Hi) can 
be represented as (|13|l for some vector ^ G Hi ® H2- In the combination 
with Lemma 3 this remark implies: 

Corollary 1. The equality (I13p establishes the isomorphism of 
Hilbert spaces Hi ® H2 and 'HS{H2, Hi). 

We now find the adjoint operator. We have: 

{^V,X) =^Tpnm{y,fnT.){en,x) = {y Jpnm {x , e„) fm) . 

nm nm 

Thus ^ 



4.2 Operation of the complex conjugation in the 
space of self-adjoint operators 

Let A G £s{W), where W is Hilbert space. We define "complex conjugate 
operator" A by its bilinear form: 

{Au,v) ^ {v,Au). (17) 

Let {fj } be an orthnormal basis in W. We find the matrix of the operator 
13 = A with respect to this basis: bij — {Afi,fj) = {fj,Afi). In the 
special case of the real basis, i.e., = /_, , we have: 

Thus its matrix is given by the transposition of the matrix of A. We 
will use the fact that the operator A is self-adjoint (we remind that A is 
self-ajoint): 

{Au,v) = {Av,u} = {u,Av) = {Av,u) = {u,Av). (18) 
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We will also use the fact that, for a positiv operator A'^, the operator A'^ 
is also positiv: 

{Ru,u) = (u^Nu) = {Nu,u). (19) 

Consider the quadratic form of the complex conjugate operator ^4 of a 
self-adjoint operator A : 

fA{4>)^CMA)-{^,M) = fA{^). (20) 

Consider the group {e, *}, where * is the operation of complex conjugation 
in a compex Hilbert space W. It induces the action in the space of real 
valued functions on : / — > /, where /(<;/>) = /(0)- We hope that the 
symbol / will not be misleading. Only real valued functions are under 
consideration. Thus it cannot be mixed with the operation of complex 
conjugation on the range of values. 

Theorem 1 implies (for A, £ Cs{H^,H,),i = 1, 2) : 

/ /Ai(<?ii)/A2(</'2)dp(0i,02) = TrDiiAi TxD22A2 + TrDi2A2D2iAi 
J H1XH2 

(21) 

4.3 The basic operator equality 

Lemma 4. Let ^' £ Hi (g) H2- Then, for any pair of operators Aj £ 
L4H,),j = 1,2, 

Tr$A2**Ai = (Ai ® A2)* = (Ai (g) A2*, *). (22) 
Proof. We have: 

Tr$A2**Ai = ^{$A2$*Aiek,ek) = X]X]V'ij(A2**Aiek,f^)(5ik 
k k ij 

= 5ZV'ij{^2**Aiei,/j) = ^^iAnjiV'i2J2(^ien.ei2)(/j2>^2/ji). 

ij «ljl '212 

By m and m we obtain (7^7, ^T") = (XT',!;;) = (fn,^/,^) = 
{A2fn,fj2)- Thus 

Tr$A2** Ai =J2J2 '/'iiji^<Aieii , ei^ ) (A2fji , fj^ ) . 

iljl i2j2 

On the other hand, we obtain: 

(Ai®l2*,*) = X]X^V»ijiV^(-4ie«i,e,,)(^2/ii,./j2) = Tr*A2**Ai. 

'131 '232 

If ^ is normalized by 1, then the right-hand side of equality ^22\ is 
nothing else than average of the observable C = Ai® A2 describing corre- 
lations between measurement of observables Ai and A2 on subsystems 5*1 
and 52 of a composite system 5* = {S\, S2) which is prepared in the state 
. On the other hand, the left-hand side of equality (|22p has the form of 
the second term in the right-hand side of formula (|10|) giving Guassian 
integral of the product of two quadratic forms corresponding to opera- 
tors A\ and A2. These mathematical coincidences provide a possibility to 
couple quantum correlations with classical Gaussian correlations, by se- 
lecting the covariance operator of the prequantum Gaussian distribution 
(corresponding to quantum state ^') in the right way. 
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4.4 Operator representation of reduced density 
operators 

Lemma 5. For any vector 4' £ H\ ® H2, the following equality holds: 

TrH2*® * = (23) 
Proof. The operator vj/v]/* acts on vector x G Hi as 
^f^'a; = ^ Vnm(a;,e„)*/m = ^^V'fciV'nm(/m,/i)(a::,e„)efc 

nm nm kl 

nk m 

Thus its bilinear form is given by 

($$*a;,y) = ^[^i)kmtpnm]{x,en){ek,y). 

nk m 

Take an orthogonormal basis {/m} in H2. The bilinear form of the 
operator TrH2 ^' ® 'l/ is given by 

(TrH2 * 8) 2/) = ^ (TrH^ * (g) *a; ® , iy « /„) 

By setting fc ^ n, Z — > fc we obtain the coincidence of two bilinear forms 
and hence the operators. 

If the vector ^l/ is normalized, then = Vl* 8 Vl* is the corresponding 



we 



density operator and pg, is the ii"i-reduced density operator. By H23I) 
obtain 

pi"=***. (24) 
f 2) 

Unfortunately, for the _ffi-reduced density operator pj^ , similar state- 
ment is not true; in general, 

pi'^ / ***. (25) 
Let find the bilinear form of the operator : H2 H2 ■ 

^*^y ^^i>nru{y,fm)^*en ^ ^ 1pnmi>ij {v , fm) {Cn , Ci) fj 

nm nm,ij 

mj i 

=^[^V™V^](y,7n:)(7^,w). (26) 

mj i 

However, for any orthonornal basis {em} in H\, we obtain: 
(TrHi* (g) = ^{-^ ®'^em®V,em®u) 

m 

= ^ (e„ ® /j ,*)(*, e„ ® fi ) {y, /, ) (/, , u) 

mjl 

/,)(/;, ^>. (27) 



11 



By setting m — J- z, j — !> m, Z — >■ j in (|27p we obtain 

(TrH,*®*y,u> =;^[;^V..V^](y,/™>(/i,^'>- (28) 

mj i 

Comparing (|26p and (|28[) we see that in general they do not coincide. To 
show this, let us take the real basis {/m}. In this case ()26[1 has the form: 

= ^[^^™V^](y,/™)(/„n>. (29) 

mj i 

Thus in this special basis the matrix elements of operators and 
TrH2*I' ® ^ are coupled via the complex conjugation. In general, we have: 

Lemma 5*. For any vector ^ G Hi (g) H2, the following equality holds 
for the operator f = Tr^/^* (g) 

f = (30) 

Proof. We now present proof which is not based on matrix elements. 
In particular, it illustrates well features of complex conjugate operators. 
For T defined in the formulation, we obtain: 



{fy,u) = {u,Ty) = {Ty,u). 

By 123 we get: 

{fy,u) = '^[^tpijtpim]{fm,y}{u, fj). 

mj i 

By ((5]) {u,v) — {v,u) and hence {u, fj) — {fj,u) and {u,v) — and 
hence {fm,y) = {y,fm)- Thus 

mj i 

In particular, if ^ is normalized (pure quantum state) then we obtain 

pi'' = (31) 

Lemma 6. Let p be a density operator in a Hilbert space W. Then, 
for any A £ £s{W), the following equality holds: 

Tr pA = Tr pA. (32) 
Proof. Let {efc} be an orthonormal basis in W. Then 

Tr pi = ^{pAck, ek) = ^(ek, pAek) 



fe fe 



k k 

Since the trace does not depend on the choice of a basis, we can select the 
real basis. Hence, it was proved that 

Tr pA = Tr Ap. 
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Finally, we will prove that 

Tr Ap = Tr pA. 

Take the basis consisting of the eigenvectors of the density operator: p = 
J2kPk^k ® Cfe. Then 

TrpA = y^(pAek,ek) = ^pk(Aek,ek) = y^(ApCk,Ck). 

k k k 

Corollary 2. Let 'i' £ Hi ® H2 be normalized (pure quantum state). 
Then 

Tr***A = Trp^^>A. (33) 

5 Classical random field description 

5.1 Ensemble of noncomposite quantum systems 

In what-follows random vectors taking values in a Hilbert space are called 
random fields. This definition is motivated by consideration of the Hilbert 
space H = L2(R.'") of square integrable functions. 

Let <7i(aj) denote a Gaussian random field in a complex Hilbert space 
H. Everywhere below we consider Gaussian random fields with probability 
distributions of the class S{H). The covariance operator of a random field 
is defined as the covariance operator of its probability distribution. 

The correspondence between QM and PCSFT in the case of a single 
quantum system with the state space H is established in the following 
way: 

1) . Density operators (statistical states of QM) are identified with 
covariance operators of prequantum random fields, p D. 

2) . Self-adjoint operators (quantum observables) are identified with 
quadratic functionals, A\-^ Ja- 

The equality ([UJ can be written as 

EfA{<t^H)^Trpl=(k),. 

It establishes the correspondence between PCSFT-averages and QM-averages. 
This story was presented in |14) - |2U| . Now we modify it. Originally the 
source of coming modification was purely mathematical - to solve the 
problem of positive definiteness in theory for composite systems, see sec- 
tion [S31 However, it happens that a natural physical interpretation can 
be provided. 

To escape measure-theoretic difficulties, at the moment we proceed 
in the finite-dimensional Hilbert space. We will come back to the real 
physical case (for infinite dimension) in section [F!3l see Proposition 4. Let 
p be a density operator. Set D — p + al, where / is the unit operator and 
a > 0. Consider the Gaussian random vector 0(a;) with the covariance 
operator D. The additional term al we can consider as (a-scaling of) the 
Gaussian normal distribution. It describes spatial white noise when the 
dimension of the space goes to infinity. 

We have EfAitj^i^)) = TrpA -I- aTrA. In this model (modification of 
PCSFT created in [15|- [20p quantum average can be obtained as a shift 
of classical average: 

{A), = EfAWco)) - aTrA. (34) 



13 



The shift is generated by the presence of the background Gaussian noise 
(say "zero point field", cf. SED, |21|-[25|). Thus QM-average can be con- 
sidered as simply normalization of average with respect to a prequantum 
random field. Normalization consists of substraction of the contribution of 
the background field. While in the finite-dimensional case the use of such 
a normalization is just a matter of test, in the infinite-dimensional case it 
becames very important. If quantum observable is represented by an op- 
erator A which is not of the trace class, then the normalization TrA = oo. 
In other words the quadratic form fA{(t>) is not integrable, cf Proposition 
4, with respect to the probability distribution pn, where p + al. Of course, 
it is a pure theoretical problem. In real experimental practice we are able 
to measure only observables represented by operators of finite ranks, see 
von Neumann |10j . Other observables (in particular, all observables given 
by operators with continuous spectra) are just mathematical idealizations. 
Nevertheless, it is convenient to have a theory which is able to operate 
with such quantities as well. From the PCSFT-viewpoint QM is such a 
theory. Thus in our approach QM has some analogy with QFT, but all 
divergences are regularized from the very beginning by choosing a special 
representation of classical averages. 

5.2 Ensemble of composite quantum systems 

Consider a composite quantum system S = (Si, S2). Here Sj has the state 
space Hj, a complex Hilbert space. Let 01 (tu) and 4>2{ijj) be two Gaussian 
random fields, in Hilbert spaces Hi and H2, respectively. Consider the 
Cartesian product of these Hilbert spaces, Hi x H2, and the vector Gaus- 
sian random field 0(aj) = (0i(aj), 02(tij)) £ Hi x H2- In the case under 
consideration its covariance operator has the block structure given by ((S)). 
Set 

(/Ai,/A2) ==-B/Ai/A2 = / /Ai(0l)/A2(02)dp(01,02). 

J HIXH2 

Set also 

COV (/Ai,/A2) = (/Ai,/A2> - (/Ai)(/A2) 

Equalities (pT|) and ((22)l imply 

Proposition 2. Let M £ Ca{Hi),i = 1,2 and let ^ £ Hi ® H2 with 
the unit norm. Then, for any Gaussian random field <j){^) Hi x H2 
with the covariance matrix D such that the non-diagonal block 

D12 = § (35) 

the following equality takes place: 

COV {fA,,fA2) = (Al ® A2*,*) = (Al ® A2)*. (36) 

This equality establishes coupling between quantum and classical cor- 
relations. In the next section we will unify classical descriptions for a 
single system, section [S] and a composite system. 

5.3 Making consistent PCSFT-models for ensem- 
bles of noncomposite and composite systems 

Operators Du are responsible for averages of functionals depending only 
on one of components of the vector random field In particular. 
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EfAii<t>i)iui)) = TrDii^i and EfA.,{(j>2){uj)) = TTD22A2. We will con- 
struct such a random field that these "marginal averages" will match those 
given by QM. For the latter, we have: 

(li)* = (Ai ® 72*, *) = Trp^''.4i, (A2>* = (Ji ® A2*, *) = Trp(,''A2, 



where 7^ denotes the unit operator in 77i,i = 1,2. By equality (|24[1 the 
first average can written as 

(li)* = Tr($$*)li. 

By equality (|32|) the second average can be written as 

and, finally, by (|31[l 

(Aa)* = Tr($**)32, 
Thus it would be natural to take 

^ _ / $$* $ 
^ 1^ $* $*§ 

Its off-diagonal block reproduces correct quantum correlations between 
systems 5*1 and 5*2 and its diagonal blocks produce correct quantum av- 
erages for system 5*1 and system 52. 

However, in general (i.e., for an arbitrary pure state *) this operator 
is not positively defined. Therefore (in general) it could not be chosen as 
the covariance operator of a random field. Let us consider a modification 
which will be positively defined and such that quantum and classical av- 
erages will be coupled by a simple rule. Thus from quantum averages one 
can easily find classical averages and vice versa. 

Proposition 3. For any normalized vector * G 77i ® 772, the operator 

( (*$* +7,/4) * \ 

""-y (*** + 72/4) J ^'^'> 

is positively defined. 

Proof. For any vector (j) — [(j>\,cj>2) G * G 77i x 772, we have: 

(L»*0,0) = ||*>l||' + J^-h(*<^2,<^l) + (**<^l,02) + ||*<?i2|P + ^l^ > 
(!|*0l||2 - ||*>l||||<^2|| + ^) + (|||<^2|P - ||'^l||||*<^2|| + ^) > 0. 

Thus operator D is positively defined|j 

We continue to proceed in the finite-dimensional case (to escape the 
problem of existence of a-additive Gaussian measure on infinite-dimensional 
space). 

For the Gaussian measure with covariance operator H37|l . we have: 
(/ai> = (^1)* - TrAi/4, (fAa) = (A2)* - TrA2/4. These relations for 
averages together with relation for correlations p6|l provide coupling be- 
tween PCSFT and QM. 

In the infinite-dimensional case Gaussian distribution with he covari- 
ance operator given by H37|) is not a-additive. 

To make it cr-additive one should consider a rigged Hilbert space: 
H+ C H C H_, where 77 = Hi x 772, and both embedding operators 
are of the Hilbert-Schmidt class. 



*Of course, the same effect can be approached by adding al for a > 1/4. 
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Proposition 4. For any normalized vector ^' G Hi ® H2, the opera- 
tor fST^I determines the a-additive Gaussian distribution on H_ or equiv- 
alently the random field <j>{ijj) valued in H_. For trace class operators 
Ai : Hi ^ Hi,i — 1,2, equalities fSd}) and flip take place. 

To prove this proposition, one should repeat the previous proofs, ex- 
istence of traces is based on the trace class condition for of operators Ai 
(and not the trace class feature of the covariance operator of a Gaussian 
measure). The crucial difference with the finite dimensional case is that 
the prequantum random field takes values not in the Cartesian product 
H = Hi O H2, but in its Hilbert-Schmidt extension. For mathemati- 
cal details, I would like to recommend the excellent short book of A. V. 
Skorohod [45| . see also |46)-|48| for applications to mathematical physics. 



6 Classical (Hilbert valued) stochastic pro- 
cess corresponding to Schrodinger's evolu- 
tion 

We again start our considerations by considering the finite-dimensional 
case. Since we do not try to go beyond QM, but only reproduce its predic- 
tions, we use Schrodinger's equation for dynamics of the "wave function'0. 
We only change the interpretation of the \I'-function of the composite sys- 
tem. Thus we start with Schrodinger's equation for a composite system 
S^{Si,S2) ■■ 

i^(t) = i?*(t), *(0) = *o, (38) 

where H is Hamiltonian of S. 

Hence, at the instant t, the covariance matrix of the prequantum ran- 
dom field (vector in the finite-dimensional case) (j){t,uj) has the form: 

jju)^ ( m)m+ii/4) *(t) 

V *W* (*(t)**(t)-f /2/4) 

The following fundamental question (having both mathematical and phys- 
ical counterparts) immediately arises: 

"Can one construct a stochastic process (valued in the Cartesian prod- 
uct Hi X H2) such that at each t £ [0, 00) its covariance matrix coincides 
with D(t)T 




6.1 Bernoulli type process 

The formal mathematical answer is yes! It is easy to construct such a 
stochastic process. Take space Q. = HtGlo 00) ^1 ^ ^2 as the space of 
random parameters, points of this space uj = {ujt) can be considered as 
functions uj : [0, cxd) Hi x H2, trajectories. Consider the family of 
Gaussian measures pt on Hi x H2 having zero mean value and covariance 
operators D(t),t G [0, oo). Consider now (on Q,) the direct product of these 
measures, -P = lite [0,00) -P*- 

Proposition 5. Let (j){t,uj) be a stochastic process having the proba- 
bility distribution P on Q,. Then, for any pair of vectors 1/1,2/2 € Hi x H2 

^So, we do not try to modify this equation, cf. |17 | 
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and any instant of time t > 0, 

E{yi,<p{t,uj)){(t>{t,u),y2) = (D(t)yi, jya), (40) 
where D{t) is given by ^4^ - 

Existence of this stochastic process is a consequence of famous Kol- 
mogorov's theorem. The equahty (|40|) is a consequence of the definition 
of probabihty P on fl. 

Thus there exists a prequantum classical stochastic process inducing 
the Schrodinger evolution for any composite system prepared initially in 
a pure state. One may say that, for a composite system, Schrodinger's 
equation describes dynamics of the nondiagonal block of the covariance 
matrix of such a prequantum stochastic process. 

This story becomes essentially more complicated after the remark that 
such a prequantum process is not uniquely determined by the D{t)\ To 
determine uniquely a Gaussian process (up to natural equivalence), one 
should define not only covariance for each instant of time, i.e, 
E{tj-i,(j)(t,u})){(f>{t,uj),y2), but so called covariance kernel D{t,s) : 

E{yi,(t>{t,ui)), (0(5,0;), 2/2) = {D{t,s)yi,y2). 

However, the formalism of QM does not provide such a possibility. It is a 
consequence of the trivial fact (but of the great importance, cf. von Neu- 
mann [10]) that Schrodinger's equation for a composite system is dynamics 
with respect to a single time parameter t, common for both subsystems, 
and not with respect to a pair of time parameters {t, s) corresponding to 
internal times of subsystems. 

Nevertheless, one may feel that the process existing due to Proposition 
5 is not adequate to the real physical situation. Since its probability dis- 
tribution P is the direct product of probabilities corresponding to different 
instances of time, it is the Bernoulli process. Its value at the instance of 
time t is totally independent from the previous behavior. Although this 
process provides right averages for each instance of time, it is hard to be- 
lieve that real physical dynamics of e.g. an electron is of the Bernoulli-type 
(and for any Hamiltonian H). We are looking for more realistic stochastic 
processes. 

6.2 Stochastic (local) dynamics in the absence of 
interaction 

We restrict our consideration to dynamics in the absence of interactions 
between and S2 after the preparation procedure. Thus we are interested 
in propagation of initially correlated random fields (vectors in the finite- 
dimensional case). Although it is a rather special dynamics, it plays an 
important role in quantum foundations. In particular, it describes the 
evolution of entanglement in the EPR-Bohm type experiments. Thus we 
consider Hamiltonian 

H ^Hi(g)l2+Ii(g)H2, (41) 

where Hj is Hamiltonian of Sj (here we use QM terminology). 
Lemma 7. Let Hamiltonian have the form J^ip . Then 

^) = e-'-^'^oe-'"^\ (42) 
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<Sf{t) = e'^^'^fo e-'"'\ (43) 

Proof. In this case 

We expand the initial state : ^'o = Tl^j i^ij^i ® fj - Then 

Thus, for 11 £ H2, we get 



By using Lemma 7 we prove: 

Lemma 8. Let the condition of Lemma 7 hold. Then 

mt)W) = e''^**;)*J,*e'^*. (44) 
^)Mt) = e'^*$J,*$oe-''^*. (45) 

Finally, we obtain: 

Lemma 9. Let the condition of Lemma 7 hold. Then the operator 
D{t) given by \39\) can be represented in the form: 



m = 



(e-'-^i'^o^o e'-"!* + J/4) 



foe 



(46) 

By using this representation it is easy to prove: 

Proposition 6. Let the condition of Lemma 7 hold. Then the operator 
D{t) given by \39\) is the covariance operator (for each instance of time 
t) of the vector process with coordinates 

Mt,^) = e-'^*eoi(^),02(t,a;) = e*^*Co2(a;), (47) 

where the initial random vector ^o{uj) = (^oi (<^), C()2(<^)) is Gaussian with 
zero mean value and the covariance operatorD{0). 

Proof. We will find not only the covariance operator for a fixed instant 
of time, but even the covariance kernel. We have, for any pair of vectors 

u,w e Hi, 

E{u,(t)i{t,uj)){(j)x{s,ijj),w) = E{e'"^^u,^oi{uj)){^oi{(^)),e'"^"w) 

= (e-'^=(ro**;, + Ji/4)e'^*u,u-). 

The same calculations can be done for the second diagonal block. Thus di- 
agonal blocks of the covariance operator of the stochastic process given by 
(|47p coincide with diagonal blocks of the operator D{t). We now consider 
nondiagonal blocks. Let now u G Hi,v G H2. We have: 

E{u,rt>i{t,Lo)){Ms,^),v) = £(e'^*u,eoi(a;))(Co2(a;)),e-'^''«) 
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The same calculations cane be done for the second nondiagonal block. 
Thus the covariance kernel has the form D{t, s) = 

(48) 

And hence, for t = s,D{t,t) = D{t). 

We emphasize that dynamics (|47p by itself is purely deterministic, 
stochasticity is generated by initial conditions. One might say that this 
process describes propagation of uncertainty of preparation. 

We remark that Propositions 5 and 6 provide two different stochastic 
processes. The covariance kernel (|48|) differs from the covariance kernel 
of the process which has been constructed by considering the product of 
Gaussian distributions pt . The latter has the covariance kernel 

D{t,s) ^ D{t)5{t - s). 

For this process, its realization at different instants of time are indepen- 
dent. The process defined by H47[l contains nontrivial dependence between 
its realizations at different times. I think that it is closer to the real 
physical situation. 

The following interesting problem arises: 

To construct a stochastic process for an arbitrary Hamiltonian, such 
that in the case of the absence of interactions this construction gives the 
process | (^7[ ). 

At the moment I am not able to solve this problem. 

6.3 Stochastic nonlocal dynamics 

We now consider another classical stochastic process reproducing dynam- 
ics of quantum correlations. 

Proposition 7. Let operator D{t) be defined by ifgffl) . Then the 
stochastic process 

at,uj) = y^VoH, (49) 
where rjQ{uj) £ Hi x H2 is distributed N{0, 1), has the covariance operator 
D{t) for any t > 0. 

Proof. Let yi , j/2 G -ffi x i/2 ■ Then 

E{yi,^{t,uj)){^{t,u;),y2) = E{y^yi,r,o{uj)){rio{u;), 

^/W)y2) = (/D(i)yi, VW)y2) = {D{t)yi,y2}. 
It is clear that the covariance kernel is given by 

D{t,s) = ^D{s)D{t). (50) 

We remind that, for Hamiltonian without interaction, we constructed 
the stochastic process <j){t,uj) given by (|47l) . In general stochastic pro- 
cesses, £,{t,uj) and (j!)(t,cj) given by (|49|) and H47p do not coincide: 

We can write the process H47|l as (j){t) — V{t)y^D{0)rjo, where V{t) = 
diag(e-*'^i,e"^2). Hence, 

D4t,s) = V{s)D{0)V{t)*. 
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On the other hand, the covariance kernel of process (|49l) is given by 
L>^(t,s) = y/D{s)D{t). We remark that D{t) = V{t)D(0)V{t)* . Hence, 
/^(I) = y(t)/^(0)V(t)*. Thus 

D^{t,s) = Vis)^/D{0)V'{s)V{t)^/D(0)V''{t). 

We remark that V*(s)V(t) ^ I,t^ s. 

The process S,{t,ijj) is nonlocal in the following sense. Its component 
£,i{t,ij) is guided not only by the Hamiltonian of 5*1, but also of 5*2; the 
same is valid for ^2{t,ijj)- Thus PCSFT (at least at the moment) cannot 
provide a definite answer on locality of the prequantum world. Quantum 
correlations can be produced by local as well as nonlocal prequantum 
stochastic processes. 

Proposition 8. In the case of Hamiltonian without interaction, see 
{7^, the stochastic process j^^P can be represented in the form: 

= e-''"'Q'l,e''^'noi + e-''"'Ql2e-''"'V02, (51) 
C2(t) = e''"'Q°2ie'"'Voi + e'"^Q°22e'''"^V02, (52) 

where 




and rio{Ld) € N{0,I). 

Proof. For example, take yi,y2 £ Hi and consider average 

S{yi,C(t,a;))(C(s,a;),2/2) 

T-i / —itH-\ itH-i / w / / \ —isH-i /-lO isH-\ \ 

= E{e "-Que '-yi,rio2(uj))(rioi(uj),e ^QiiC ^^2) 
+E{e'''^-[Ql2re'"'V^,V^2{^)){m2{i^),e'"^{Ql2re"^'y2) 
+E{e-""^ QW»^ yi , r?io(^)) (r?2o(^) , e""- (Qja)*^^^'^^ 2/2) 
+E(e''"^- {Ql2re'"' y, , 7710 (a;)>(,72o(a;), e^^^ Ql^e^'"' ^2) . 
Two last terms are equal to zero, since 

E{zx,r[ia{bj)){r]2o{^),Z2) = 
for any pair z\ £ H\, 22 £ H2- The first two give us 

/ —itH-\ y^O itH-i —isH-i ^0 isH-\ \ 

(e iQiie ^2/1,6 ^Qne ^^2) 

I / itHo / \* itH-i isHo / r^O \* isH-\ \ 

+ {e ^{Qi2) e ^1/1,6 ^(Qi2) e 17/2). 

Thus 

7-1 /i \ —isHt y-\0 i(s — t)H-i /^O itH-i , —isH-i ^0 i(t — s)H'y //— \0 \* itH-] 

Dii{t,s) = e ^Qiie^ ' Que (5i2e ^ ' -(Qi2)e ^ 

Representation (|51ll . (|52|l implies that even in the absence of interaction 

between the subsystems Si and S2 of the system 5* the dynamics of 5*1 
depends on the Hamiltonian H2 and vice versa. It can be interpreted as 
a sign of "action at the distance". Thus the stochastic process ^(t) can 
be considered as "nonlocal" - opposite to the process ^(t) given by (|47|) . 
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6.4 Infinite-dimensional case 



To proceed in the infinite-dimensional case, one should consider a rigged 
Hilbert space: H+ C H C H_, where H — Hi x H2, and both embedding 
operators are of the Hilbert-Schmidt class. Stochastic processes take val- 
ues in the Hilbert space _ff_ (and not H — Hi x H2)- All previous results 
are valid for any unitary dynamics ^'(t) = U{t)'^o. 

This paper was written under support of the grant "Mathematical 
Modeling" of Vaxjo university and the grant QBIC of Tokyo University 
of Science. It was presented at "Feynman Festival", June, 2009; the au- 
thor would like to thank Vladimir Manko for his critical comments which 
improved understanding of the model. 
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